Modification of smeared phase transitions by 
spatial disorder correlations 



David Nozadze, Christopher Svoboda, Fawaz Hrahsheh, and Thomas Vojta 

Department of Physics, Missouri University of Science and Technology, Rolla, MO 65409, USA 

Abstract. Phase transitions in disordered systems can be smeared if rare spatial regions develop 
true static order while the bulk system is in the disordered phase. Here, we study the effects 
of spatial disorder correlations on such smeared phase transitions. The behaviors of observables 
are determined within optimal fluctuation theory. We show that even short-range correlations can 
qualitatively modify smeared phase transitions. For positive correlations (like impurity atoms attract 
each other), the order parameter is enhanced, while it is suppressed for repulsive correlations 
(like atoms repel each other). We use computer simulations to generate various types of disorder 
correlations, and to verify our theoretical predictions. 
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INTRODUCTION 

Phase transitions occur between different macroscopic states of many-particle systems 
when control parameters such as temperature, pressure, magnetic field, chemical com- 
position, etc. are varied. Phase transitions occurring at nonzero temperatures are called 
thermal or classical transitions. They are driven by thermal fluctuations. Another type 
of transitions, the so-called quantum phase transitions occur at absolute zero tempera- 
ture, when a nonthermal external parameter is changed [1]. They are driven by quantum 
rather than thermal fluctuations. The study of observable quantities at quantum phase 
transitions is an important subject of modem condensed matter theory. 

Since materials often feature considerable amounts of quenched disorder, the study 
of phase transitions in the presence of imperfections or randomness has received much 
attention. Initially, it was thought that any kind of disorder can destroy sharp phase 
transitions because (in the presence of disorder) the system divides itself up into regions 
which can undergo the transition at different values of the control parameter Later, 
it became clear that phase transitions are generically sharp in the presence of weak 
disorder. However, under special conditions, rare finite-size spatial regions can indeed 
undergo the phase transition independently, leading to a smearing of the transition. This 
happens, for example, at quantum phase transitions with overdamped order parameter 
dynamics [2, 3, 4] or at thermal phase transitions with extended defects [5, 6]. 

Theories of disordered systems often assume the randomness to be uncorrected 
in space even though real systems contain some amount of disorder correlations. At 
critical points, this assumption is justified as long as the correlations are not too long- 
ranged. This can be understood, e.g., by looking at the Harris criterion [7] for the 
stability of a clean critical point against weak disorder. Harris divided up the system into 



regions whose linear size is the correlation length . The variation of the local critical 
temperatures from region to region can be estimated using central limit theorem, yielding 
5Tc ^ ^ '^Z^, where d is the space dimension. Harris observed that the critical behavior 
of the clean transition can be unchanged if the fluctuations dTc between the regions 
are smaller than the global distance t = T — Tc from criticality, i.e., dT^ < t. Because 
t, ~ this leads to the inequality t^^l^ < t for ? — > 0, where v is the correlation length 
critical exponent of the clean system. Thus, for stability of the clean critical behavior we 
must have dv > 2. This inequality is called the Harris criterion. 

Weinrib et ah, [8] generalized the Harris criterion to long-range correlated disorder. 
They considered power-law correlations , C(x — y) = [(71.(x) — Tc){Tc{y) — 71.)]av ~ |x — 
y|^", where [...]av is the disorder average. To estimate fluctuations between correlation 
volumes i^^, the central limit theorem estimate has to be replaced by 

^^Tcf=^ I, I, Vc(x - y) ~ — "(«'^) . (1) 

For a — d, we find {5Tc)^ ~ ^^log^. Then, the stability criterion for clean critical 
behavior takes the form - 2 > for a > J or av - 2 > for a < J. 

Thus, for a >d, the original short-range Harris criterion is recovered. However, for 
a <d, long-range correlations lead to a new inequality. Thus, as long as the correlations 
decay faster than x^'^ with distance, correlated and uncorrelated disorder have the same 
effect on the stability of the clean critical behavior. 

In Ref. [9], we studied the effects of spatially correlated disorder on smeared phase 
transitions, and we showed that even short-range correlations qualitatively modify the 
behavior of observable quantities. Here, we summarize this theory. We then report the 
results of new computer simulations for which we have generated spatially correlated 
atom distributions using Kawasaki Monte Carlo simulations [10]. 

CORRELATED DISORDER AT SMEARED PHASE TRANSITIONS 

In this section, we summarize the theory of Ref. [9] for the behavior of observables in 
the tail of smeared phase transitions. We consider a material consisting of substances A 
and B. Material A is in the magnetic phase characterized by a negative distance from 
criticality ?a < 0. Material B is in the paramagnetic phase with ?b > 0. The quantum 
phase transition from a non-magnetic to a magnetic phase can be tuned by substituting 
magnetic atoms A for nonmagnetic atoms B in the binary alloy Ai^^Ex- 

Due to statistical fluctuations in the atom distribution, there can be large (rare) spatial 
regions which show local order because they have a high local concentration of A 
atoms even if the bulk system is in the non-magnetic phase. If the order parameter 
fluctuations are overdamped, the quantum dynamics of these regions freezes. This leads 
to inhomogeneous magnetic order and a smearing of the phase transition [2]. 

We now explore the effects of disorder correlations on smeared phase transitions 
within optimal fluctuation theory. Assume that there are attractive short-range correla- 
tions between like atoms, with disorder correlation length (^dis. Like atoms tend to form 
clusters of typical volume Vdis ~ 1 +^"3^8 • Roughly, the transition point where the 



phase transition would happen without having any rare regions effects, can be found by 
setting the average distance from criticality to zero. This gives jc^ = — ?a/ (^b — ^a)- A rare 
region of linear size Lrr can show local static order if the local concentration of non- 
magnetic B atoms is smaller than some critical concentration Xc (Lrr) which can be esti- 
mated from finite size scaling [ 1 1 ] , jc^ = -'J^c — ^^Lrr ■ Here, D is a non-universal constant, 
and is the shift exponent which takes the value of 1/2 within mean-field approxima- 
tion. As the concentration Xc must be positive, we obtain the condition Lmin = (Z)/jc^) ^/"^ 
for the minimum size of a locally ordered rare region. 

The probability P{N, Nb) for a rare region with N = L^^ sites to have Nb = A'^xioc sites 
occupied by B atoms is equal to the probability i'cius(A'^/Vdis, A^s/Vdis) of finding A^s/Vdis 
B-clusters among N /V^i^ clusters. It is given by the binomial distribution 

^clus(A^/Vdis,^loc/Kus) = L /y^ J (1 ''"^ • (2) 

To calculate the total magnetization, we can simply integrate the binomial distribution 
(2) over all rare regions showing local order. In the regime where the concentration x is 
somewhat larger than x*^, we find (up to power-law prefactors) 



M ~ exp 



c (x-xQ)2-^/^ 



(3) 



where C is a non-universal constant. In this regime, the magnetization drops exponen- 
tially with X, and spatial disorder correlations modify the prefactor in the exponent. 

In the far tail of the smeared transition at jc ^ 1, the composition dependence of 
the order parameter is given by a non-universal power law M ~ {1 —x)^ with j8 = 
(^^min + ^^^s) /^dis- For Small disorder correlation length ^^{^ <ti Ln^n, earlier results 
for uncorrelated disorder [12] are recovered, /3 ~ L^j^- In the limit of large correlation 
lengths ^ ^ L^in, all clusters show local order, and all magnetic A atoms contribute 
to the magnetization. This leads to j8 = 1. Thus, in the far tail, disorder correlations 
modify the exponent of the order parameter. In both regimes, for attractive correlations 
as considered here, the magnetization M at a given concentration x increases with 
increasing correlation length ^dis- 



COMPUTER SIMULATIONS 

To illustrate the theory we now show results of numerical simulations of a toy model. 
Motivated by the quantum-classical mapping [1] we consider a classical Ising model 
with infinite-range interaction in one (timelike) dimension and short-range interaction 
in three spacelike dimensions. The Hamiltonian has the form [9] 

1 1 

H=-— ^ /o'S'y.T'S'z t' - — ^ /ySy^TrSy -;:/ . (4) 
(y,z),T,T' '^y,T,T' 

Here, y and z are space coordinates, T is the timelike coordinate. is the system size 
in the time direction, (y, z) denotes pairs of nearest neighbors in space. Jy is a binary- 



random variable, which takes values i/ or Jj^ depending on the type of atom on the given 
lattice site y. The concentration of sites with 7/ and are x and l—x, respectively. The 
values of Jy at different lattice sites are not independent, but they are correlated. 

Because the interaction is infinite-ranged in the timelike direction, this dimension can 
be treated within mean-field theory. This leads to coupled mean-field equations for the 
local magnetization at site y 

my = tanh j8 {Jyniy -\- Lz■^'^z + h) , (5) 

where h= 10~^^ is a small symmetry-breaking field and the sum is over all nearest 
neighbors of site y. P = l/Td is the inverse classical temperature (not related to the 

physical temperature of the quantum system which is encoded in L^) [1]. The smeared 
phase transition can be tuned by changing the composition x in the temperature range 
Th > Tci > Ti, where Tfi—Jh + 6Jo and Ti — J[ + 6Jq are the phase transition temperatures 
for pure systems with all Jy = Jf, and Jy=Ji, respectively. 

In order to generate correlated binary random variables we u se a more realistic method 
than the one used in Ref. [9]. We consider a model of A and B atoms with short-range 
interactions. This model is equivalent to an Ising model 

H' = -Jkw "y"z > (6) 

where ny = ±1 corresponds to A and B atoms and Jkw > for attractive correlations. 
To obtain correlated atom distributions with a well-defined concentration x, we use 
Monte Carlo simulations by the Kawasaki method [10]. We can control the correlations 
between the atoms by changing the "Kawasaki" temperature T^w in these simulations. 
At Tkw = 4.5, atoms are correlated with power-law correlation, while at higher Tkw 
correlations are short-ranged (decay exponentially). Correlation length increases from 
to °o with decreasing Tkw from oo to 4.5. Figure 1 show examples of the atom 
configurations for two Tkw and several concentrations x. 

To study the smeared phase transition we then solve the mean-field equations (5) nu- 
merically. Figure 1 show the magnetization M as a function of concentration x for several 
disorder correlations (several Tkw)- The data show that at a given x, the magnetization 
M increases with decreasing Tkw from oo to 4.5 (i.e. with increasing correlation length). 
Thus, the magnetization is enhanced compared to uncorrelated case. We also consider re- 
pulsive (anti-)correlations. To obtain repulsively correlated atoms, we set Jkw < when 
generating the atom distribution. The solution of (5) shown in Figure 1 demonstrates 
that the magnetization is suppressed compared to the uncorrelated case. 

CONCLUSIONS 

We have studied the effects of spatial correlations between impurity atoms on smeared 
phase transitions. We have shown that even short-range correlations between atoms have 
dramatic effects and qualitatively modify the behavior of observables. Attractive disor- 
der correlations enhance the magnetization in the tail of the smeared transition, while 
repulsive correlations suppress it. This is in contrast to conventional critical behavior, at 



FIGURE 1. Left: Examples of the atom configuration on a square lattice of 96^ sites generated by 
Kawasaki Monte Carlo simulations with attractive interactions Jkw = 1 ■ The Kawasaki temperature Tkw 
is 4.5 in the top row and 32.0 in the bottom row. The concentration x takes values 0. 1, 0.3 and 0.5 (left to 
right). Right: Magnetization M versus composition x for correlated disorder ranging from strong attractive 
correlations to strong repulsive (anti-)correlations. The data stem from one disorder realization of 256^ 
sites for each case, using 7/, = 20, Ji = %, Jo = I, h = 10^'", and T^-i = 24.25. 

which correlations that decay sufficiently fast lead to the same critical behavior as un- 
correlated disorder. We have verified our theoretical predictions by performing computer 
simulations. To this end, we have generated both attractively and repulsively correlated 
disorder by means of Kawasaki Monte Carlo simulations. We have found that the re- 
sulting data qualitatively agree with our theory. We emphasize that our results do apply 
not only to the quantum phase transitions considered here, but to all disorder induced 
smeared phase transitions. 
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